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ABSTRACT 

At high redshift, the universe was so young that core-collapse supernovae (SNe) are 
suspected to be the dominant source of dust production. However, some observations 
indicate that the dust production by SNe is an inefficient process ^ casting doubts on 



the existence of abundant SNe-dust in the early universe. Recently, iPerlev et al.l ([20101) 
reported that the afterglow of GRB 071025 - an unusually red GRB at z ~ 5 - shows 
evidence for the SNe-produced dust. Since this is perhaps the only high redshift GRB 
exhibiting compelling evidence for SNe-dust but the result could easily be affected by 
small systematics in photometry, we re-examined the extinction properties of GRB 
071025 using our own optical/near-infrared data at a different epoch. In addition, we 
tested SNe-dust models with different progenitor masses and dust destruction efficiencies 
to constrain the dust formation mechanisms. By searching for the best-fit model of the 
afterglow spectral energy distribution, we confirm the previous claim that the dust in 
GRB 071025 is most likely to originate from SNe. We also find that the SNe-dust model 
of 13 or 25M0 without dust destruction fits the extinction property of GRB 071025 
best, while pair- instability SNe (PISNe) models with a 17OM0 progenitor poorly fit 
the data. Our results indicate that, at least in some systems at high redshift, SNe 
with intermediate masses within 10 — SOMq were the main contributors for the dust 
enrichment, and the dust destruction effect due to reverse shock was negligible. 

Subject headings: gamma ray : burst — dust, extinction — galaxies: high-redshift — 
galaxies: ISM 



1. Introduction 

Dust plays an important role in understanding the galaxy formation and evolution over the 
cosmic history. Dust absorbs and scatters ultraviolet (UV) and optical lights from young stars. 
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making it difficult to comprehend how stars were formed in galaxies using the optical and the near- 
infrared (NIR) observations that sample the redshifted UV/optical lights of high redshift objects. 
Although the infrared and submm observations allow us to get a glimpse of the light from young 
stars re-emitted in the rest-frame infrared (e.g. Choi et al. 2006), such observations for the high 
redshift universe are still technically challenging. Therefore, it is imperative to know how the dust 
extinguishes the UV/optical light at high redshift, in order to understand how stars and galaxies 
were formed in the early universe. 

In the local and the intermediate re dshift (z < 3) universe, A symptotic Giant Branch (AGB) 
stars are one of the major dust sources (jFerrarotti &: Gail 1120061 ). When the universe was much 
younger (~ 1 Gyr) or at z > 5, however, AGB stars cannot be the main producers of dust simply 
because stars at that epoch do not have enough time to evolve into the AGB-phase. 



I t is also known that core-collapse supernova e (SNe) can produce dust grains (jKozasa et al 



19891 : iTodini &: Ferraral l200ll : iNozawa et al.ll2003l ). Since core-collapse SNe can occur on short 



time scales, the dust in the early universe is speculated to be formed from SNe (IMaiolino et al 



2004 ) . However, the dust production from c ore collapse SNe is known to be inefficient (j Green et al 



2004 : lErcolano et al.l 120071 : iRho et al.l l2009l ) , making it unclear if core-collapse SNe can produce a 
sufficient amount of dust to significantly obscure star formation activities at high redshift. 

Observations of broad absorption line quasars at high redshift and some gamma ray bursts 
(GRBs) suggest that their rest-frame spectral energy distributions (SEDs) can be fitted best if 
one assumes a dust extinction curve produced by SNe-dust (Maiolino et al. 2004; Hirashita et al. 
2005, hereafter H05; Stratta et al. 2007; Gallerani et al. 2010). In such a study, GRBs can be 
greatly useful tools. A GRB afterglow spectrum in UV/optical follows a simple power law that 
results from the synchrotron radiation of accelerated relativistic electrons by internal and external 
shocks propagating into the surrounding medium of GRB progenitors, making it easy to trace dust 
extinction curves imprinted on the SEDs . The extreme brightness of GRBs helps us obtain the 
SEDs of them at high redshift even with a small telescope. 



Currently, observational constraints on dust in high redshift GRBs are weak. While IStratta et al, 

(120071 ) claimed that the extinction prope rty of GRB 05090 4 at z = 6.29 is due to SNe-dust, 
a re-analysis of the GRB 050904 data by I Z afar et al.l (|2010l ) reveals no strong evidence to sup- 
port the SNe-dust. Other high redshif t GRBs are also known to harbor little amount of dust 



( Cucchiara et al 



2011 



Zafar et al 



2011 



)• 



GRB 071025 is an unusually red GRB believed to be at a redshift of z ~ 5 (Im et al. 2007; 
Perley et al. 2010, hereafter PIO). Its re d color and the inf lexed shape of the afterglow SED suggest 
dust extinction dominated by SNe-dust (jPerlev et ahlboiol l. Other than GRB 050904, GRB 071025 
is currently the only GRB which is claimed to harbor SNe-dust. However, as we saw from the case 
of GRB 050904, subtle systematic errors in photometry could lead to a large difference in the 
derived dust extinction properties (Zafar et al. 2010; PIO). Therefore, in order to provide an 
independent assessment of the dust extinction properties of GRB 071025, we analyzed the SED 
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of its afterglow using our own dataset of optical/NIR images taken in RIJHKg and constructed 
at a different epoch from PIO. In addition to the verification for SNe-d ust, we also examined a 
variety of SNe-dust extinction models from H05 and lHirashita et alj (|2008l ) (hereafter H08) to place 
constraints on the progenitor mass and the importance of the destruction of the dust by external 
medium. 



2. Analysis 

2.1. Observation and Data Analysis 

We started the follow-up imaging observation of the afterglow of GRB 071025 using the 1- 
m telescope at Mt. Lemmon Optical Astronomical Obse rvatory (LOAO) on 2007 October 25, 
04:26:54.07 UT, ~1080 sec after the Swift BAT trigger (jPagani et al.ll2007l l. We obt ained 



V, R and I band im aging data till the end of October 25, and a part of October 26 (jim et al 



20071 : iLee et al.ll20ld ). Our data showed detections of a bright afterglow in the / band and a faint 
counterpart in the R band. The afterglow was not detected in the B and V bands. We reduced 
the data using the standard reduction procedures of dark subtraction and flat-fielding using IRAF 
packages. We obtained the photometry of the afterglow and surrounding stars using SExtr actor 
with the parameters of 5 p ixels for DETECT MIN A REA, 1.2 a for DETECT THRESH, and 1.2 a 
for ANALYSIS THERSH (iBertin &: Arnoutall996l ). The magnitudes are computed using circular 
apertures with an aperture radius of 1.5 FWHM of point spread functions and aperture corrections 
derived from bright stars around the afterglow which amount to 0.2 - 0.4 mag. We obtained the 
data for standard stars (SA95 and SA115) in the second night which were used to do the photometry 
cahbration of the GRB 071025 field. 

The J, H, and Kg band data were obtained by the Simultaneous Quad Infrared Imaging Device 
(SQIID) on the Kitt Peak May all 4-meter telesco pe at Steward Observatory, with the start time of 
2007 October 25, 05:56:59 UT (j.Iiang et ahlboOTl ). The total exposure time in each band was 240 
second, and the three bands data were taken at the same time. The photometry was calibrated 
using 2MASS point sources brighter than J = 16.0, H = 15.5, and Kg = 14.0 mags within 120 
arcsec distance from the afterglow. Additional J, H, and Kg data were obtained at At ~ 1.05 days 
with the Canada-France-Hawaii Telescope. 

The results of the photometry are presented in Table 1, where the magnitudes are given in the 
AB system. For the conversion of the V ega-based magnitu des into AB magnit udes, we used the 
conversion relation of lFrei &: GunnI (|l994l ) for RI bands, and lCiliegi et al. I (bood) for JHKg bands. 
Galactic extinction is also corrected using E{B — V) = 0.07 ( Schlegel et al.l 1998 ). which amounts 
to 0.19, 0.14, 0.06, 0.04, 0.03 mags in RIJHKg, respectively. 
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2.2. Synchronization of Photometry 

The afterglow data were taken at different epochs but GRB afterglows change their brightness 
and in some cases their colors, making it necessary to construct a SED at a given epoch for the 
subsequent analysis of the afterglow. For this purpose, we use the mid-time of the JHKg data (6605 
sees after the burst trigger) as the epoch for which we construct the SED, since the JHKg data 
were taken simul taneously. Note that we use external shock model with the form of Ft^v ~ 



(jSari et al.lll998l ) for the GRB 071025 afterglow, where t is the delay time since the BAT trigger. 



1/ is the obseved frequency, a is temproal index and (3 is spectral index. 

Before determining the SED at this epoch, we consider if it is appropriate to consider such a 
SED as a simple power law SED from synchrotron radiation of afterglow by analyzing their temporal 
evolutions. At t > 1500 sec, PIO find that the light curves in the X-ray, and the optical/NIR are 
described by simple power-laws. B y fitting single po wer-law functions to the light curves of the 



publicly available Swift XRT datarl ijEvans et al.ll2007l ) and to the I, J, and Ks-hand data including 



those presented in PIO, we find that over an interval oft = 1900 to 17000 sec ax = 1.75±0.06, a/ = 
1.50 ±0.05, aj = 1.44 ±0.04, and or = 1.41 ±0.03, or aopUcai/NiR- ax ^ -0.25 to -0.34 ±0.08. 
When the analysis is extended to i ~ 10^ sec, we get aj — ax = —0.27 ± .04. These result s, 



being consistent with aopUcai — cnx < —0.25 of uniform ambient environment (jUrata et al.l 120071 ) . 
suggest that cooling break (vc) lies between the X-ray and the optical wavelengths and that the 
light in optical/NIR is dominated by the afterglow from external shock. Since we also find no late 
time prompt signal (e.g., flare and shallow decay) in X-ray, we conclude that stable afterglow light 
characterizes the SED at t = 6605 sec. 

The remaining task is to synchronize R and / band photometry to this epoch. One of the I 
band images from LOAO is nearly contemporaneous to the JHKg data, lagging behind the NIR 
observation only by 209 sec. To correct for the small difference in the observed time, we interpolate 
the I band data with the power law model. 

We fit the / band light curves only with our own data and with additional / band data from 
PIO. The fitted values of a are 1.48±0.11 and 1.48±0.04 without or with the PIO data respectively. 
The necessary correction to derive the contemporaneous / band magnitude from the 6814 sec post- 
burst data point is ±0.050 mag with a = 1.48. Because the choices of a in both two cases share 
the same value in spite of different errors, we will adopt the correction of 0.05 mag to compute the 
/ band magnitude at 6605 sec, which gives / = 19.18 ±0.16 mag. For R band photometry, we used 
two methods. First, we extrapolated the R band detection at t = 1500 sec tot = 6605 sec assuming 
an achromatic decay of the light curve during this period as found in PIO. This gives the R band 
photometry of R = 21.3 mag. As the second approach, we derived a limit in the R band photometry 
from a stacked image made of seven 300 sec R band frames around t = 6605 sec. In the stacked 
image, no afterglow was detected, giving an upper limit in R band photometry at 21.1 mag at 3 a. 



^http: //www. swift. ac.uk/xrt_curves/00295301/ 
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These two methods give consistent results, therefore we adopt R = 21.3 mag for the R band pho- 
tometry at t = 6605 sec. For the 1 a uncertainties of data, 0.17 mag is used in / band to reflect the 
uncertainty in a. In the case of R band, we adopt the error, 0.3 mag, of the R frame at t = 1500 sec. 



3. Analysis Method 



The observed photometry data are fitted using model SEDs as described below. The under- 



lying afterglow emission is assumed to follow a simple power law of F^, 



as produced by 



synchrotron emission of relativistically mo ving par t icles. The attenuation of photons below Ly— a 
is implemented following the prescription of lMadaul (|l995l ). For the intrinsic dust ext i nctio n of GRB 
071025, we try the M ilky Way (hereafte r MW ), LMC, SMC extinction laws (jPeil Il992l ). and an 
extinction curve from ICalzetti &i KinnevI (|1994| ) which is usually applied to star- forming galaxies. 
For SNe-dust extinct ion curves, we tried an empirical extinction curve derived from z ~ 6 quasars 
(jMaiolino et al.ll2004l ). and theoretical extinction curves of H05 and H08 (see discussion in Section 
4.2 for more detail of the model). The m odel SEDs were convoluted with the transmission curves 
of filters and quantum efficiency of CCD (iLee et al.ll2010l ). 



4. Resutls & Discussion 

4.1. Redshift of GRB 071025 

We derive the redshift of GRB 071025 before proceeding to a detailed comparison of the best- 
fit models using different extinction curves. By fixing the redshift, the direct comparison between 
different models becomes more straightforward. Fixing the redshift to a single value can be justified 
if the derived redshift is not strongly dependent on models. 

In order to determine the redshift, we fit the observed SED with different extinction curves 
with four free parameters: redshift. Ay, beta, and the normalization factor. We find that the 
derived redshifts are all betwen 4.6 to 4.85 due to the sharp break at R band with R — I ^ 2.0 
which can be naturally explained with intergalactic absorption of the light below Ly-a. The result 
is consistent with the value reported in PIO. 

If we change the redshift 2 a apart from the best fit result for each extinction curve, the 
amounts of increased x^/^ values for SN-dust extinction curves approximately range from 3 to 5 
: 3.48 (Maiolino), 4.21 (Hirashita 25M0), and 5.46 (Hirashita 13Mq) whereas other extinctions 
show higher values than SN-dust extinctions do : 7.90 (SMC), 9.12 (LMC), and 12.40 (Calzetti). 
The MW extinction shows a larger increase, 16.23 compared to the others, which can be explained 
by the infiuence from 2175A bump. Within 2 a deviation of redshit from the best fit, SN-dust 
curves always present less x^/'^ than the other extinctions. Consequently, none of all extinctions 
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have better fitting results when the redshift deviated by 2 cr from the best fitting one is chosen. 
Since the dependence of the redshift on the assumed dust extinction curves is small, we will fix 
the redshift of GRB 071025 to be a fiducial value of 4.8 in the following analysis of the extinction 
property. We also tried the same analysis using different values for the redshift, and found that the 
conclusion reached in the paper is not affected by the choice of the redshift. The difference in the 
reduced values between different models does not 



4.2. Dust Properties 

Figure 1 shows the result of the SED fit with various dust extinction curves (SMC, MW, 
Calzetti, and Maiolino). Like in PIO, we observe an inflextion in the observed SED shape - As 
we go from the shorter wavelength to the longer wavelength, IJHKg data points move up and 
down with respect to a simple power law model SED. While the best-fit models with SNe-dust 
have the \^ /v values less than 1, models with the dust extinction curves of LMC, SMC, MW, 
and the Calzetti produce the \^ /v values of about 4 — 5. The MW and the LMC extinctions are 
characterized by a 2 175 A bump which is absent in the SED of GRB and the best-fit Ay^s of these 
models converge to zero. The SED with the Calzetti curve is featureless, and similar to the SED 
with no extinction. Although the adoption of the SMC extinction curve improves the goodness of 
the fit a little bit compared to the MW or the Calzetti curves, it cannot reproduce the observed 
changes of slopes between JHKg bands as well as the SNe-dust extinction curves. 

Since some of the fits converge to no extinction, we tried a fit without dust extinction, which 
has a benefit of reducing one free parameter. The fit returns the of 2.55, an improvement 

over the fits with the MW and the LMC extinction curves. However, the best-fit spectral slope 
is too steep at /3 = 1.6 in comparsion with fix ^ 1 observed by the XRT, which is inconsistent 
with the synchrotron theory, (iopticai,iR = /?x — 0.5 if Uc is between X-ray and optical/IR. As our 
light curve analysis suggests that the afterglow emission follows the external shock model with 



3(p-l) p-1 . 

F{t,v) ~ t 4 1/2 for the optical/NIR where p is electron power law index (jSari et al.Ml998l ) 



From our optical/NIR light curves, we get p ~ 2.87 and consequently /3 ~ 0.94. Also, GRB 



afterglows are known t o have (3 ~ 0.6 and /3 < 1.1 for all the well-studied cases (jKann et al.l 120061 . 



2010l : IZafar et ahlboilh . Therefore, we argue that no extinction model is unfavorable 



The SNe-dust models perform far better than the models with the other dust extinction laws, 
confirming the previous result from PIO. 



4.3. Progenitor and Environment of the Dust-Producing SNe 

The H05 and H08 allow us to test various SNe-dust production scenarios with a broad progen- 
tior mass range from 13 to 17OM0, with or without mixing of heavy elements inside the He-shell, 
and with or without dust destruction due to the reverse shock. Dust grains can be destroyed by the 
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SNe shock, modifying the grain size distribution and eventually the extinction curve. In H08, the 
shock destruction of the dust within SNe is implemented, and the destruction efficiency depends on 
the ambient gas density and it is varied from 0, 0.1, and 1 cm~^, with being the case for no dust 
destruction and 1 cm~^ for the most efficient dust destruction. If the gas density is much larger 
than 1 cm~^, the dust destruction is so efficient that there would be almost no dust extinction 



(jNozawa et al.l 120071 ) . Figure 2 and Table 2 summarize our finding after the GRB 071025 data 



fitted with these models. 

First, we discuss which of the mixed or the unmixed core model fits the data better using 
25Mq progenitor models as references. The right panel of Figure 2 shows that the unmixed core 
model fits the data better than the mixed core model. Mixed core only produces oxidized molecules, 
since most carbons are locked in CO. Extinction is dominated by Si02 in such cases, regardless of 
progenitor mass. Since the Si02 ~ dominated extinction curve monotonically increases with 1/A 
without any change of the slope, it cannot explain the infiexing SED shape of GRB 071025. We 
find that the same argument applies to models with different progenitor masses, therefore, our 
discussion to unmixed core progenitors in the following. 

Next, we examine if the dust destruction is needed to account for the observed SED. The 
triple dot-dashed and the dotted lines represent models with dust destruction at ni^=0.1 cm^'^ and 
1 cm~^ respectively. The destruction affects small grains, and the effect is reflected in the model 
SED as shown in the left panel of Figure 2. After the destruction, the extinction curve becomes 
featureless resembling the SMC or the Calzetti curves, and the afterglow model SEDs with the 
dust destruction effect fit the data as poor as those with the SMC or the Calzetti curve. Therefore, 
we conclude that dust destruction should be negligible for GRB 071025. We checked that this 
conclusion holds for models with other progenitor masses. 

Finally, we present what the most likely progenitor mass is for the SNe responsible for the dust. 
For this, we compare the best-fit models with 13, 20, 25, 30, and 170Mq progenitors with unmixed 
core and no dust destruction (right panel of Figure 2). Figure 2 shows that the models with the 
progenitor mass of 13 or 25M0 fit the data significantly better than the others (the 'X^ jv < 1 
versus ~ 5). This is because the infiexing SED shape between / and Kg bands can be explained 
by extra absortions at the rest-frame 1200 — 2000A and at around 3000A , which results from 
the increased contribution of carbon grains (13Mq model), or Mg2Si04 and FeS grains (25M0 
model) with respect to Si grains. The 30Mq progenitor model fails to fit the data because of the 
excessive increase of Si grains. The 17OM0 model, given as a represen tative case for Pair Insta bility 



Supernovae (PISNe) which could be abundant in the early universe (jHeger &: Wooslevi |2002| ) . also 
gives poor fits to our data. This results from the high contribution of Si grains that make the curve 
look more or less monotonically increasing, and the deficiency of the extra extinction from Fe that 
introduces a plateau in the model SED at shorter wavelengths. 

Additionally, in order to confine the most optimal progenitor mass of SN-dust extinction for 
GRB 071025, we introduce the Y band data from PIO. We interpolated the Y data from REM 
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and MAGNUM in PIO to determine Y magnitude at t = 6605 sec through ?. Although the Y 
magnitude was inserted into the SED, both Hirashita 13 and 25Mq models show acceptable fitting, 
as x^/^ = 2.48/(5-3) and 0.57/(5-3) for 13 and 25Mq, respectively. When the dust destruction is 
taken into consideration, the Y magnitude do not still enhance the poor fits of models where dust 
destructions exist. In conclusion, the additional information of Y do not alter the consequences of 
our own dataset for appropriate SN-dust extinctions. 

5. Summary 

We examined the dust property of a red afterglow of GRB 071025. Analysis of our own dataset 
which is largely independent from PIO supports the evidence for the existence of SNe-dust at z ~ 
5. Furthermore, using the SNe-dust extinction curves of H05 and H08, we find that models with 
13 and 25Mq progenitors with unmixed core and no dust destruction give the best fit to the data. 
Given uncertainties in the model, we suggest that the most plausible progenitors of the SNe-dust 
are intermediate mass stars with 10 — SOMq. 
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Fig. 1. — SEDs of GRB 071025 at the 6605 sec post-burst epoch, fitted with models incorporating 
various extinction curves. The soUd hnes indicate the best-fit model curves, while the dashed line 
shows a result with no extinction. The Maiolino model, given as a representative of the SNe-dust 
extinction curves, fits the data best. The dotted lines for each panel are plotted with no extinction 
assumption. 
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Fig. 2. — Comparison of the GRB 071025 SED with SNe-dust models with different physical 
ingredients such as dust destruction, progenitor mass, and mixed/unmixed core. In the left panel, 
we show 25M0 progenitor models with mixed core (the long-dashed line), with unmixed core (the 
solid line), and with unmixed core with dust destructions where the ambient densities are nff=0.1 
cm~^ and nH=l-0 cm~^. In the right panel, we indicate cases of different progenitor masses with 
unmixed cores under no dust destruction. In both figures, the dotted lines indicate cases without 
dust extinction which are virtually identical to cases with very efficient dust destructions. The figure 
shows that the models with unmixed cores, 13 and 25M0 progenitors, and no dust destruction fit 
the data best. Y magnitude is extracted from PIO. 



Table 1. Optical & NIR Photometry of GRB 071025 



Filter Facility Mid-point time elapsed Total integration time Magnitude "'^ 
after Swift alert 





LOAO 


1091 s 


30 s 


17.33 ± 0.08 




LOAO 


1901 s 


300 s 


17.27 ± 0.09 




LOAO 


3703 s 


300 s 


18.46 ± 0.09 




LOAO 


5061 s 


300 s 


18.81 ± 0.15 




LOAO 


6605 s 




19.18 ± 0.16 




LOAO 


6814 s 


300 s 


19.23 ± 0.16 


R 


LOAO 


1563 s 


300 s 


18.95 ± 0.17 


R 


LOAO 


6605 s 




21.3 ± 0.30 


J 


SQIID 


6605 s 


240 s 


18.18 ± 0.05 


H 


SQIID 


6605 s 


240 s 


17.83 ± 0.06 


K 


SQIID 


6605 s 


240 s 


17.31 ± 0.05 


J 


CFHT 


90589 s 


1200 s 


22.13 ± 0.35 


H 


CFHT 


92308 s 


1050 s 


21.56 ± 0.24 


K 


CFHT 


88792 s 


1200 s 


20.78 ± 0.12 



^The magnitudes are in AB magnitude system 
^Galactic extinction was corrected 
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Table 2. Fitting Results Varying Extinction Curves 



Extinction curves " 




^3000 




None 


1.56±0.11 




5.10/(4-2) 


Milky Way 


1.56±0.11 


~0 


5.10/(4-3) 


LMC 


1.56±0.11 


~0 


5.10/(4-3) 


SMC 


1.05±0.47 


0.30±0.24 


3.87/(4-3) 


Calzctti 


~0 


2.57±0.07 


3.94/(4-3) 


Maiolino SNe 


1.24±0.17 


0.77±0.34 


0.17/(4-3) 


Hirashita 13M0(U) 


0.86±0.32 


1.73±0.77 


2.48/(5-3) 


Hirashita 20 & 3OM0(U) 


1.54±0.10 


~0 


7.29/(5-3) 


Hirashita 2OM0(M) 


1.19±0.37 


0.15±0.15 


6.44/(5-3) 


Hirashita 25M0(U) 


0.79±0.30 


2.65±1.00 


0.57/(5-3) 


Hirashita 25M0 (M) 


1.26±0.23 


0.07±0.05 


5.40/(5-3) 


Hirashita 25M0 nH=0.1{\J) 


1.36±0.20 


0.77±0.73 


6.37/(5-3) 


Hirashita 25M0 n//=1.0(U) 


1.53±0.16 


0.09±0.88 


7.26/(5-3) 


Hirashita 17OM0(U k M) 


1.54±0.10 


~0 


7.29/(5-3) 


Hirashita ITOM© nif=0.1(U) 


1.54±0.10 


~0 


7.29/(5-3) 


Hirashita ITOM© nif=1.0(U) 


1.54±0.10 


~0 


7.29/(5-3) 



^U:Unmixed core, M:Mixed core 



